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HERMITIAN SYMMETRIC SPACES OF TUBE TYPE 
AND MULTIVARIATE MEIXNER-POLLACZEK POLYNOMIALS 

Jacques Faraut & Masato Wakayama 



ABSTRACT. Harmonic analysis on Hermitian symmetric spaces of tube type is a natural 
framework for introducing multivariate Meixner-Pollaczek polynomials. Their main properties are 
established in this setting: generating and determinantal formulae, difference equations. As an 
^sO ', application we consider the problem of evaluating moments related to a multivariate Barnes type 

integral involving the Harish-Chandra c-function of a symmetric cone. 

< 

u 

^ ' In this paper we revisit the identity 

a' 



(i±|) n = i + 2 c(m, n)x m+ \ 



m=0 



where 

OS 
(N 



c(m ' n)= £C)G+i 



for n E Z. The coefficients c(m, n) appear in the formula giving the moments of the 
statistical distribution \i n of the eigenvalues of a random Hermitian matrix in case of the 
Gauss unitary ensemble (GUE): 



m m (»n) = - \ tr(X 2m )F n (dX) : 

n JHerm(n,C) 



where P n is a Gaussian probability on the space Herm(n, C) of n x n Hermitan matrices: 

F n (dX) =C n e- tr{x2) X(dX) 
(A is the Lebesgue measure). In fact 

1 (2m)! 
9tt2m(>n) = -^z — rc(m, n). 
n 2 m m! 

See [Harer-Zagier,1986], §4, [Mehta,1991], formula 5.5.30, [Haagerup-Thorbj0rnsen,2OO3], 
§4, [Faraut, 2004], §5. We will see how this identity appears in the framework of harmonic 

1 



analysis on the unit disc and upper halfplane, and that the coefficients c(m, n) are essen- 
tially Meixner-Pollaczek polynomials. In [Davidson-Olafsson-Zhang,2003] a generalization 
of this formula appears in the framework of Hermitian symmetric spaces of tube type. We 
will consider such an extension by using a different method. 

Let us describe the scheme of the paper in the one variable case. We start from 
weighted Bergman spaces on the unit disc. The Cayley transform maps these spaces to 
weighted Bergman spaces on the right half plane. Then, by an inverse Laplace transform, 
we obtain an L 2 -space on the positive half-line. Finally we perform a Mellin transform and 
get an L 2 -space on R. Simultaneously, starting from the monomials on the disc, we obtain 
on each step an orthogonal basis, and, at the last step, the Meixner-Pollaczek polynomials 
as, essentially, Mellin transforms of Laguerre polynomials. In this way, starting from the 
binomial formula, seen as a generating formula for the monomials, one obtains a generating 
formula for the Meixner-Pollaczek polynomials. Similarly, from the fact that the monomials 
are eigenfunctions of the Euler operator, it follows that the Meixner-Pollaczek polynomials 
are eigenfunctions of a difference operator. This scheme extends almost word for word in 
the framework of Hermitian symmetric spaces of tube type, and leads to a natural definition 
of the multivariate Meixner-Pollaczek polynomials and their properties. In this extension 
the positive half-line becomes a symmetric cone, and the Mellin transform generalizes as 
a spherical Fourier transform. 

In the first section we consider the one variable case, to make clear the scheme which 
is carried on in the general case. In Section 2 we develop our construction in the setting 
of Hermitian symmetric spaces of tube type. The invariant differential operators on a 
symmetric cone turn out to be a powerful tool. Multivariate Laguerre polynomials, which 
appear in the third step of the construction, have been considered with the same point of 
view in [Davidson-Olafsson,2003], [Aristidou-Davidson-Olafsson,2006] and [2007] (see also 
[Chebli-Faraut,2004]). Eventually multivariate Meixner-Pollaczek show up as, essentially, 
spherical Fourier transforms of multivariate Laguerre polynomials. Similar analysis has 
been considered in [0rsted-Zhang,1994], related to invariant differential operators on a 
Hermitian symmetric space. This construction leads to a generating formula, as it is shown 
in Section 3. In case of the multiplicity d = 2, we establish in Section 4 determinantal 
formulae for multivariate Laguerre and Meixner-Pollaczek polynomials. We show in Section 
5 that multivariate Meixner-Pollaczek polynomials are solutions of a difference equation, 
and satisfy a Pieri's formula. In the last section an application is given to the evaluation 
of multivariate moments related to a Barnes type integral involving the Harish-Chandra 
c-function of a symmetric cone, yielding an alternative proof of a result by Mimachi [1999a]. 

It is quite natural to develop this analysis in the framework of the special functions 
associated to root systems. This point of view is considered in [Sahi-Zhang,2007]. 

Part of this work has been done as the first author was visiting the Faculty of Mathe- 
matics of Kyushu University. This author would like to thank the colleagues of this faculty 
for their hospitality. 

Mathematics Subject Index: 32M15, 33C45, 43A90 

Keywords: Meixner-Pollaczek polynomial, Laguerre polynomial, Hermitian symmetric 
space, Jordan algebra, spherical function. 
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1. Functions spaces on the unit disc and the complex half-plane, one vari- 
able Meixner-Pollaczek polynomials. — In the first part of the paper we will consider 
the one dimensional case, and develop some classical analysis on weighted Bergman spaces 
on the unit disc V in C, and on the right halfplane T, and then about Laplace integral 
representation of holomorphic functions and Mellin transform. Finally Meixner-Pollaczek 
polynomials will show up in this setting. 

(1) For v > 1 let Hl(V) denote the Hilbert space of holomorphic functions / on the 
unit disc V = {w G C | \w\ < 1} such that 

11/112 = 4 1} ! l/H! 2 (l - \w\ 2 r- 2 m(dw) < oo, 

where m is the Lebesgue measure on C. The constant a v = is choosen such that the 
function <po = 1 has norm 1. The monomials <p m (w) = w rn form an orthogonal basis of 
Hl(V), and 

1 1 i ll2 "7,! 

\y )m 



(2) The Cayley transform 

1 + w 



w i— > z = c(w) = 



1 — w 



maps the unit disc V onto the right half-plane T = {z = x + iyeC\x>0}, and its 
inverse is given by 

-if \ z ~ 1 
c (z) = . 

V ; z + 1 

For a holomorphic function / on £>, define the function F = C v f on T by 



'w-^w-( £ T i n(Si) 



The map C v is a unitary isomorphism from H 2 {V) onto the space H 2 {T) of holomorphic 
functions F on T such that 



\F\\l = a {2) [ \F{x + iy)\ 2 x v - 2 dxdy 
Jt 



< oo. 



The constant ai ' is such that the function Fq , Fq'(z) = ( £ ^-) , image of <f>o, has 
norm 1. The functions 



form an orthogonal basis of 7^(T), and 



\ F (»)\ 



ml 



(3) The functions in H^(T) admit Laplace integral representations. Define the modi- 
fied Laplace transform C v of a function ip on ]0, oof by 



/•OO 

Jo 



with ai, 3 '* = j^y- For G 1^(0, oo), equipped with the norm given by 



/» OO 

^0 



the function £ u ifj is holomorphic on T, and £^ is a unitary isomorphism from L^(0, oo) 
onto TiliT). In particular, if tpo(u) = e~ u , then £ !y 'i/ , o = Fq, and ?/>o nas norm 1. If 
ip(u) = e~ u u k , then 



v+k 



nv poo / 2 \ 

(A,V0(*) = fH i, e " (z+1) ^ I/+fc " 1 ^ = 2_fe ^ fe (iTI ) 

Define the Laguerre functions as 

^ ) («) = e-^- 1) (2n), 

where Lm denotes the classical Laguerre polynomial of degree m: 

7>-i)M = ^ r(m + z/) (-x) fc = H»^/mW , 
" + fc!(m-fc)! m! gU/M* 

Hence 

#(M)= (^g ( _ 2)k 0J_ e -^. 

Let us compute the modified Laplace transform of the Laguerre functions: 

fc=0 v 7 

by the binomial formula. It follows that the Laguerre functions ipm form an orthogonal 
basis of L„(0, oo), and that 

ll«H 2 ,= (//) " ! 



m! 
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(4) Finally we perform a Mellin transform. We define the modified Mellin transform 
of a function ip on ]0, oof as 

1 f°° 

MMs) = -, u \ / ^{u)u s ^- l du. 
By the classical Plancherel theorem 

J mu)\\^du=-J |A^zA)| 2 |r(zA+-J| 2 dA, 

and M. v is a unitary isomorphism from L 2 (0, oo) onto 

L 2 (M, M(dA), with M(dA) = a^ } \T(iX+^)\ 2 d\. The constant a[, 4) is such that the function 
M-u^q = 1 has norm 1: 



4" - 1-*- 



2tt T(u) ' 

If ip(u) = e~ u u k , then .M^^) = (s + § )fc is a polynomial of degree fc. We define 
the polynomials qm as the modified Mellin transform of the Laguerre functions tf)m^'- 
q^} = M. u ipm . The polynomials qm , when restricted to the line iR, are orthogonal 
with respect to the measure M(d\). Hence their zeros are located on the imaginary axis. 
Furthermore 

W\\l'-= r \q m {iX)\ 2 M{d\)= { ^. 

J — oo rn - 

The polynomials qm are essentially Meixner-Pollaczek polynomials. Recall the hy- 
pergeometric representation of the Meixner-Pollaczek polynomials (see [Andrews- Askey- 
Roy,1999], p.348): 

P«(A; <P) = ^hle im * 2 F 1 (-m, a + zA; 2a; 1 - e" 2 ^). 

Proposition 1.1. 

q^(s) = ^ 2 F 1 (-m,s+^u;2) 1 or ? &>(iA) = {-i) m P m (A; |) . 
Furthermore one can write 



Proof. Recall that 



k=0 



m + u-l\ fs- § 
m — k J \ k 
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Therefore 

tf(s) ^| ( _ 2)i ( : )^ (s+i) ^Mp 2Fl( _ m , s+i; , ;2) . 

Observe further that 

(•+5)* = (-D*« ("V S )- □ 
This construction provides a generating formula for the Meixner-Pollaczek polynomials 

q y m ■ 

Proposition 1.2. — For z/, s e C, iu e C with \w\ < 1, 

00 

E ?ff(«v m = (i - + «o— *• 

m=0 

Proof. We may assume Re z/ > 0, and Re(s + ^) > 0. The statement for u, s G C will then 
follows by analytic continuation. Define, for \w\ < 1, \(\ < 1, 

m=0 

By the classical binomial expansion, 

H( 1 \w,0 = (l-wC)-". 
We apply now the transform C u to both handsides with respect to ( and get 

^WhE^I^K", ( 2 ) 

m=0 

and 

Hj?\w,z) = 2"(1 - w)- v {z + c(w))~". 
Next we apply the inverse of the modified Laplace transform C v with respect ot z: 

00 

H < ?\w,u)=Y J ^\u)w m , (3) 

and 

Hj?\w,u) = {l-w)- u e- u ^. 
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Finally we apply the modified Mellin transform M. v with respect ot u: 

oo 

H^(w,s)=J2tiHs)w™, (4) 



m=0 

and 



H?\w, s) = (1 - w) s -% (1 + «;)" s -^ D 

Since Hl 4 \—w, s) = Hc, 4 \w, —s), it follows that has the parity of m: qm\—s) = 
( l) m ^ } (s). 

For z/ = this gives 

1 oo 

(i^) = £ 2 fi(-m, a; o;2V», 



m=0 

and, with s = —n, one gets the identity in the introduction: 



LXJ 



m=0 



In this setting we obtain also a difference equation for the Meixner-Pollaczek polyno- 
mials q^} , and a three terms relation. Let D v be the following difference operator: 

D v f{s) = {s+ U -)f{s + l)-{s- U -)f{s-l). 

Proposition 1.3. - The Meixner-Pollaczek polynomials q$ are eigenf unctions of 
the operator D u : 

D v q%> = (2m + v)q£>, 
and satisfie the following three terms relation: for m > 1, 

2sq^\s) = (m + v - \)q^_ x {s) - (m + l)^*). 

Proof, a) Let us consider on the disc D the modified Euler operator Dv = 2w-^ + v. 
The monomials (j) m (w) = w m are eigenfunctions of it: Dl 1} (/) m = (2m + v)(f> m . 

Then on the half-plane T consider the operator such that C v D\p = DvC v . 
Then 

D?) = (z * _ + vz, D^F^ = (2m + „)i#\ 



since = C u (p r 
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Further consider the operator Djp on ]0, oof such that C v = C U D^\ Then 
= -u^ + = (2m + v)^\ 

since C v i>m = -Fm ) • 

Finally let the operator D { „ ] on C be such that M v Df ] = M v . Then = D u , 
and, since M. u ipm = q™, , we get 

D v q£> = {2m + v)q£>. 
b) Proposition 1.2 can be written, with s = m + ^, 

oo 
£;=0 

Let us apply the operator to both sides. Since D^\w k ) = (k + u)w k+1 — kw k ~ 1 , and 
Di 2) F{n } = (2m + v)Fm \ we get, for k > 1: 

(2m + z/)^ } (m + = (fc + - 1)^ (m + - (k + 1)^ (m + . Q 

Notice that the operators D$ (i = 1, 2, 3, 4) are essentially selfadjoint in corresponding 
Hilbert spaces (with appropriate domains). 

The fact that the zeros of the polynomials qm are located on the imaginary axis is 
related to the so-called Local Riemann hypothesis (see [Bump-Choi-Kurlberg-Vaaler,2000]). 
In a recent paper by Kuznetsov ([2008]) the three term relation is used to study the 
expansion of the Riemann S-function in series of Meixner-Pollaczek polynomials, and the 
zeros of the partial sums are investigated. 

2. Function spaces on a Hermitian symmetric space of tube type and 
multivariate Meixner-Pollaczek polynomials. Let V be an irreducible Hermitian 
symmetric domain of tube type realized as the unit ball in a simple complex Jordan algebra 
Vc of dimension iV and rank n, seen as the complexification of a Euclidean Jordan algebra 
V. The Cayley transform 



w i — ^ z — c(w) = (e + w)(e — w) 



-l 



maps the domain V onto the tube domain Tq = {z = x + iy G Vc | x G O}, where O 
denotes the symmetric cone in V. Let G be the identity component in the group G(Q) of 
linear automorphisms of O, and K C G the isotropy subgroup of the unit element e G V . 
The space V(Vc) of holomorphic polynomials on Vc decomposes multiplicity free as 

V{V c ) = @V m 

m 



where V m is a finite dimensional subspace, irreducible under G. The parameter m is a 
partition m = (mi, . . . , m n ) G N n , mi > • • • > m n . The subspace of i^-invariant poly- 
nomials in V m is one dimensional, generated by the spherical polynomial $ m , normalized 
by the condition <fr m (e) = f . We consider on V(Vc) the Fischer inner product, and denote 
by K m (wi, W2) the reproducing kernel of V m - (See [Far aut-Koranyi, 1994], Chapter XL) 

The Gindikin gamma function Tn of the symmetric cone O is defined, for s G C n , 
with Resj > |(j - 1), by 



r n (s)= / e- tr ^A s (x)A(x)-^m(d, 
Jo. 



'x), 

where A is the determinant on Vc, A s is the power function defined on and the multi- 
plicity d is such that N = n + ^n(n — 1). (Notation m will denote the Euclidean measure 
on any Euclidean vector space.) Its evaluation gives 



n j 

N — n -1— r _ , CI 



i=i 

hence Tq extends analytically as a meromorphic function on C n . 

(See [Far aut-Koranyi, 1994], Chapter VII.) For a G C, and a partition m, one writes 
a + m = (a + mi, . . . , a + m n ) G C, and defines the generalized Pochhammer symbol 

r^(q + m) 

r^(a) 

For v > 2f - 1, Hlip) denotes the Hilbert space of holomorphic functions / on T> 
such that 



I = a U) / |/(«;)| 2 /i(«;)-- 2 v m (^) < oo. 
Jv 



Recall that h(w\, W2) denotes a polynomial on Vc x Vc, holomorphic in w±, antiholomorphic 
in W2, such that, for u G V, h(u,u) = A(e — -u 2 ), and, for w G Vc, h(w) = h(w,w) by 
definition. The constant ai} is such that the function $0 = 1 nas norm 1. The spherical 
polynomials $ m form an orthogonal basis of the space 7i 2 {T>) K of i^-invariant functions 
in Hl(D), and 

1 (-) 

ll/f) ||2 V n ) m 

dm (^)m 

where <i m = dim"P m . 

For a function / holomorphic on V, one defines the function F = C u f on Th by 



F(z) = {C v f){z) = A( Z -±^y"f((z -e)(z + e)" 1 ). 

The map C v is a unitary isomorphism from Hl(V) onto the space 7i 2 (Tn) of holomorphic 
functions F on Th such that 



F|| 2 = [ \F(z)\ 2 A(x)"- 2J ^m(dz) < 00. 
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The constant affl is such that the function = C u $o, Fq V \z) = A( £ ^ £ ) has norm 

1. The functions Fm' = C v Q> m form an orthogonal basis of HI(Tq) k of if -invariant 
functions in H^(Tq). 

The functions in 7Y^(Tq) admit a Laplace integral representation. The modified 
Laplace transform C v given, for a function ip on fi, by 

= «i 3) / e- (z|u V(«)A(w) i/ -»m(^), 
in 

is an isometric isomorphism from the space L^(fi) of measurable functions ip on O such 
that 

|2 



2 = a (?) /* |^( w )| 2 A( w )"-f m (^) < 00, 
Jo, 



onto H^iTct). The constant is such that the function ^o(u) = e _tr ( w ) has norm 1, and 
then £ u ^o = ^b- Let us recall the definition of the multivariate Laguerre polynomials. 
The generalized binomial coefficients (™) are defined by the generalized binomial formula: 



kCm ^ ' 



(Observe that (™), as a function of m, is a shifted Jack polynomial. With the notation of 
[Okunkov-01shanski,1997], 

m^_P k *(m;0) Q= d 



see formula (3.8).) We define the multivariate Laguerre polynomial Lm -1 ' by 



^(->-^E(-)^>- 

V / m 



(This definition slightly differs from the one given in [Far aut-Koranyi, 1994], p. 343), but 
it has the advantage that it extends exactly the one variable case.) Define the Laguerre 
functions as ^\u) = e" tr ^ L^ _1) (2u). Then 

V n I m 

([Faraut-Koranyi,1994], Proposition XV.4.2.) The Laguerre functions form an orthogonal 
basis of Ll(Q) K of if -invariant functions in 1^(0). 

Finally we perform a spherical Fourier transform which can be seen as a multivariate 
Mellin transform. Recall that the spherical function <p s , for s G C n , is defined on O by 



<p s (u) = / A s+p (ku)dk, 
Jk 
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where p = (pi, . . . , p n ), pj = j(2j — n — 1). (dfc denotes the normalized Haar measure on 
the compact group K.) We introduce now the modified spherical Fourier transform T v as 
follows: for a function ip on f2, 

J- n(s + 2 J 

Observe that J- v ^q = 1. By the spherical Plancherel theorem, T v is an isometric iso- 
morphism from L^{Vt) K onto the space L 2 (R n ,M(d\)) " of square integrable symmetric 
functions on IR n with respect to the measure 

M(d\) = \T Q (i\ + p+ ^) | 2 ^J_m(dA), 

where c is the Harish-Chandra c- function and ai^ is such that M(d\) has total measure 
1: 

! \F u iP{i\)\ 2 M{d\) = a { ^ [ \ifj(u)\ 2 A(u) u -^m(du). 
The symbof of a partial differential operator D on V is defined by 

DeW® = a D {x, 0e (x|?) (x, £ E V). 

(D acts on the variable x.) If D e D(O), the algebra of G-invariant differential operators 
on ~ G/K, then <7£) is a G-invariant polynomial on V x 1/ in the following sense: for 
g E G, aD(gx,£) = (Td(x, <7*£). The map D \— > o~d(x, e) is a vector space isomorphism 
from D(O) onto the space V(V) K of K-invariant polynomials on V. The spherical function 
cp s is an eigenfunction of every D e D(O): D</? s = 7r>(s)<^ s , and the map D i— > 725 is an 
algebra isomorphism from D(O) onto the algebra V(C n ) &n of symmetric polynomials on 
C n . 

In particular, let D m denote the invariant differential operator D for which ctd(x, e) = 
$ m (x), and let 7 m = 7D m - Observe that, by the spherical Taylor formula ([Faraut- 
Koranyi,1994], p.244), 

-Tk(m-p). 



k ^ (?)„ 

Lemma 2.1. - Letip(u) = e~ tr ^p(u) , where p is a K -invariant polynomial on V. 
There is D e D(O) swc/i £/iai ctd(x, — e) = TTien V'(' u ) = -De _tr<u ' ) , and £/ie modified 

spherical transform of tp is given by 

F u ip(s) = 7D (-s- |). 
Hence q = T v i\) is a symmetric polynomial on C n . 
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Proof. First part of the statement follows from what was explained above. Let us perform 
an integration by part as follows: 



/ e~ u{u) p(u)tp s (u)A(u)%~^m(du) 
= [ D(e- tr ^)ip s+ »(u)A(u)-^m(du) 
= [ e- tr(tt) (DV s+! )(«)A(«)-^m(cl«) 



= 7^*(s+^) / e- tr M<p s+ »(u)A(u)-%m(du) 

= 7o-(s+|)r n (s+|), 

and recall that 7d*(s) = 7r>(— s) ([Faraut-Koranyi,1994], Proposition XIV. 1.8). ;1 

We define the multivariate Meixner-Pollaczek polynomial Qm (s) as the modified 
spherical Fourier transform of the Laguerre function ^m" 1 : 

Qfc ) (s) = Jv*S ) (s). 

Observing that, by Lemma 2.1, 

^(e- tr ^ k )(s) = (-l) |k| 7k(-s-^), 

we obtain: 

Theorem 2.2. - The multivariate Meixner-Pollaczek polynomials Qm form an 
orthogonal basis of the space L 2 (R n , M(dA)) ™. Furthermore 

U/m kCm 7 

Remark. In the one variable case (V = R), 7fc(s) = [s]fc, and 7fe(— s — f) = (— l) fc (s + 

3. Generating formulae. - Following the same scheme as in Section 1, we will 
obtain the following generating formula for the multivariate Meixner-Pollaczek polynomi- 
als: 

Theorem 3.1. - For w e V, v e C, s e C n ; 

^mQ^sj^H = A(e - w 2 )-%<p s ((e - «;)(e + w)- 1 ). (4) 

m 

As in the one variable case, it follows that Qm\— s) = (— l)' m 'Qm^(s). In [Davidson- 
Olafsson-Zhang,2003], the polynomials p^m are defined through this generating formula. 
The definition on p. 179 is slightly different: 

4Q^(s) =Pi/,m(is). 
12 



Proof. The reproducing kernel K v for the weighted Bergman space Hl(D) is given by 

K v (w t w') = h(w,w')- u = J2(u)mK m {w,w'). 

m 

We consider its average over the compact group K, letting w' = (: 

H^(wX)= [ h(wM)- v dk. 

JK 

It is holomorphic in w and £. Then we obtain 



m V n ) m 



(1) 



We apply the transform C v to both handsides with respect to (. The left handside 
becomes 

Hi 2 \w,z) =2 n "A(e -w)~ v [ A(kz + c(w))~" dk. 

Jk 

We have used the identity A(z + e)h(w,c~ 1 (z)) = A(e — w)A(z + c(w)). Since CV$ m = 
Fm\ we obtain 

HW(w,z) = J2d m -^^ m (w)F^(z). (2) 

m V n / m 

Next we apply the inverse of the modified Laplace transform L v with respect to z. 
Define, for w £ D, u £ fi, 



\c(w)) 



H^(w,u) = A{e-w)~ u ( e~( ku ^ w, )dk. 

Jk 

The modified Laplace transform of Hv(w,u) with respect to u is equal to h1 2 \w, z). In 
fact 

<ynv r 

___ J e -^Hl 3 \w,u)A(u) u -^m(du) 
2 



n 



^—A(e-w)- u [ ([ e-^ k ^e-( ulc(w) )dk)A(uy-%m(du) 
2 nu A{e-w)- v [ A(kz + c{w))~ u dk = H^\w,z). 

JK 



On the other hand, since: 

M 



n ) m 

we get the generating formula for the Laguerre functions: 
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Proposition 3.2. 

^^ffM*-^) =H<?\w,u). (3) 



Finally we perform the modified spherical Fourier transform T v with respect to u. 
Define, for w G D, s G C r , 

i/W («,, s ) = A(e - «; 2 )-^^ s ((e - w)(e + w)' 1 ) 
= A(e-w 2 )-%p- s (c(w)), 

since </? s (;r _1 ) = </?_ s (x). The modified spherical Fourier transform of hI s \w,u) with 
respect to u is equal to Hl 4 \w, s). In fact 

l n(s + j J ■'fi 

= J_-^A(e -«;)-" / e-(" l(e+u,)(e -" rl )v? s (w)A(w)^-^m(^) 

= A(e - w)"> s+f ((e - iu)(e + w)" 1 ) 
= A(e - w 2 )~^ s ((e - w)(e + w)' 1 ). 

By definition the Meixner-Pollaczek polynomial Qm is the modified spherical Fourier 
transform of the Laguerre function . Hence: 

J2dmQ { ZHs)$ m (w) = A(e - w 2 )-^ s {(e -w)(e + w)- 1 ). (4) 

m 

D 

4. Determinantal formulae. - In the case d = 2 (V = Herm(n, C), K = U(n)), 
there are determinantal formulae for the Laguerre functions ^m"* and for the Meixner- 
Pollaczek polynomials Qm . Consider a Jordan frame {ci,...,c n } in V, and let 5 = 
(n-l,n-2,...,l,0). 

Theorem 4.1. - Assume d = 2. The multivariate Laguerre function admits 
the following determinantal formula involving the one variable Laguerre functions : 

forU = YTj = l u i c i> 

m ~ V(ui,...,u n ) 

where V denote the Vandermonde polynomial: 

V(lii, . . . ,U n ) = - Mj)- 

14 



As a result one obtains the following determinants formula for the multivariate Laguerre 
polynomials: 

K(til,...,U n ) 

9/. We start from the generating functii 

3.2) 



Proof. We start from the generating function for the Laguerre functions \l/m^ (Proposition 



Hl 3 \u t w)=A(e-w)- v [ e -(ku\(e+wKe-w) ^ 
m 

In the case d = 2, the evaluation of the following integral is classical: for x = Yl7=i XiCi > 
V = T^=iVjCj, then 

Z(x,y) J^e dk 5. _ _ _ ^ x n )V(y 1 , . . . , y n ) ' 

In fact X{x,y) is the so called Itzykson-Zuber integral: 

l(x,y) = [ e tT{uxu * y) du. 

JU(n) 

Therefore, for u = Ya=i u i c i> w = Y!j=i w i c h 

detfe Ull ~ w i ) 

HW(u,w) = S\T\(l-w j )- u ^-ttx tx — 7- 

j=l V ^ Ul ' •••i u "n)V \ i- Wl > • • • > 1- Wn ) 



Noticing that 



1 + Wj 1 + W k _ Wj - w k 



1-Wj l~W k (1 + Wj)(l+W k )' 

we obtain 

_ 1 + w j \ 

1 det((l-«; 7 )" (!/ " n+1) e J 

(«, «,) = 5!2 - 2 ^ 

V[ui,...,u n )V{w 1 ,...,w n ) 

We will expand the above expression in Schur function series by using a formula due 
to Hua. 

Lemma 4.2. — Consider n power series 

00 

fi(w) = J2 c( S wm (< = l,-..,n). 

m=0 
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Then 

det(fi(wj)) ^ 
V{w 1 ,...,w n ) ^ 

where s m is the Schur function associated to the partition m, and 

a m = det(cf£. +s .). 

(See [Hua,1963], Theorem 1.2.1, p.22). 

Let v' = v — n + 1, and consider the n power series 



00 



fi(w) = (1 - w)- v 'e- u ^ = ^\u i )w m . 



m=0 



Since 



we obtain 



n 

dm^m^^WjCj^j =S m (w 1 ,...,W n ), 



3 = 1 



gfffr) = m-hjn-l) ^^ ("*)) Q 

V(ui,...,u n ) 

By using the same method we will obtain a determinantal formula for the Meixner- 
Pollaczek polynomials Qm . 

Theorem 4.3. 

det(q^ n + + s 1 ( Si )) 
V ^Si, . . . , s n j 

where qm denote the one variable Meixner-Pollaczek polynomials. 

Proof. We start from the generating formula for the Meixner-Pollaczek polynomials Qm 
(Theorem 3.1): 

A(e - w 2 )~^ s ((e -w)(e + w)- 1 ) = ]T d m Q^ (s)^ m (w). 

m 

For x = Ym=i x i°i, the spherical function (p s (x) is essentially a Schur function in the 
variables xi, . . . , x n : 



(p s {x) = 5l(xix 2 ■ ■ .x r ) 2 



I(n-l). det ^) 



V(s 1 ,...,s n )V(x 1 ,...,x n )' 
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Let us compute now, for w = Y^=i w j c ji 



A(e-w 2 ) *(p s ((e-w)(e + w) x ) 



11 llM + WiJ v( Sl ,..., Sn) v(^,...,^) 



3 = 1 3 = 1 ~ ' 1 



In the same way as for the proof of Theorem 4.1, we obtain 



A(e-w 2 ) %(p s ((e-w)(e + w) x ) 
= 5!(-2)-' 



^n(n — 1) V /_ 



V(si,...,s n )V(«;i,...,ii; n ) 
We apply once more Lemma 4.2 to the n power series 

oo 

fi(w) = (1 - u;)«-t(1 + = ^C( S ,K 

m 

with v' = v — n + 1, and obtain finally: 



<&'(,)^(-2)-t^ y^ □ 

K(Si,...,S n ) 



5. Difference equation and Pieri's formula for the multivariate Meixner- 
Pollaczek polynomials. - - Let us first recall Pieri's formula for the spherical functions. 



Proposition 5.1. 



(tTx)ip s (x) = ^2aj(s)p s+Sj (x), 

3 = 1 



Sj Sfc 



with 

,w = n 

({sj} denote the canonical basis ofC n .) 

See [Dib, 1990], Proposition 6.1 (with a minor correction), where it is called Kushner's 
formula. (See also [Zhang, 1995], Theorem 1). One observes that 

c(s) 



(Xj{s) = 



c(s + e 3 y 
17 



in agreement with the asymptotic behaviour of the spherical function <p a \ 

for s = (si, . . . , s n ) real, with s\ < . . . < s n , and a = Yl7=i UiCi w ith a\ < . . . < a r , 

V? s (exp to) ~ c(s)e (p+s|a)t (£ -> oo). 
Recall that the c-function is a product of beta functions: 

tf [Sj - Sk, 

j<k 



c(s) = c Jl b(sj -s k , 



By puting m = s + p (recall that pj = |(2j — n — 1)), Sj = — |j + |(n + 1), one 
obtains 

n 

(trx)$ m (x) = ^Oj(m)$ m+ej .(:r), 
i=i 

a (m) = TT mj -^"f^-^" 1 ) 
" A " rrij - m k - - k) 



with 



(In agreement with Lassalle's results [1998], p. 320, 1.-4.) 

We introduce the difference operator acting on functions on C n : 

(A,/)(s) = + \--^ n ~ I))",' (s)/(s + e,) 

i=i 

n v d 

3=1 

Theorem 5.2. - (i) The Meixner-P ollaczek polynomials Qm are eigenf unctions 
of the operator D v : 

A,QM = (2|m| + i/)Q&>. 
(ii) Furthermore they satisfy the following Pieri's formula: 

n 

^2(mj + u-l--(j- l))aj(m - - p)d m - ej Q^_ e . (s) 



" 5Z( m ^ + 1 + ^ ~ ^')) a i(- m ~ £ J + P) rf m+ ej <5mi ej (S). 

J = l 

Statement (i) has been obtained in another way in [Davidson,6lafsson-Zhang,2003]: 
Theorem 6.1. 
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Proof. (1) We consider first the differential operator on the bounded domain V as 

D^f(w) = 2(w, V/H) + nuf(w). (1) 

This is a modified Euler operator. (We denote by (z, w) the C-bilinear form on Vc extending 
the Euclidean inner product on V so that (z\w) = (z,w).) The functions $ m , being 
homogeneous, are eigenfunctions of D„ 

D^<S> m = (2|m| + nv)$ m . 

(2) Next we define the operator on the tube domain Th as the image of 
through the Cayley transform. More precisely D^C V = C V D^}\ 

Proposition 5.3. 

D^F(z) = (z 2 - e, V/(z)) + utr(z)F(z). (2) 

(2) (2) 

This means that the symbol o v of D„ is equal to 

al 2 \zX) = (z 2 -e,() + vtrz. 



Proof. We will use the following formulae: V(A(x) a ) = aA(x) ot x~ 1 (Proposition III. 4. 2 
in [Faraut-Koranyi,1994]), and, if j denotes the inversion: j(x) = a; -1 , then (Dj) x = 
— P{x~ x ) (P is the so-called quadratic representation of the Jordan algebra V). By writing 
c(w) — 2(e — to) -1 — e, we obtain the differential of the Cayley transform: (Dc) w = 
2P((e - w)- 1 ). If F = C v f, i.e. f(w) = A(e - w)~ u F (c(w)) , then 

DPf(w) = 2(Vf(w),w) + ™>fW 

= 2z/((e - w)~\ w)A(e - w)- u F(c(w)) 

+ 4A(e - w)-"(VF(c(w)),P((e - + nuf(w). 

With z = c(w), one gets 

2((e - w)- 1 , w) = tr(z) - n, 4P((e - = - e. 

Therefore 

D^f(w) = A(e-w)- v G(c(w)), 

with 

= (z 2 -e,VF(z)) + vti(z)F(z). D 



The functions Fm are eigenfunctions of : 

D^F^ = (2\m\ +nu)F^. 
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(3) Now the differential operator Dl is defined on O as the inverse image of D v 
under the modified Laplace transform. 

Proposition 5.4. 

(3) 

This means that Dl is the differential operator with symbol 
Proof. By writing 

= («- ( ' h "P) iS(0) . 

and observing that is formally selfadjoint, we obtain 

c(D^)w = (^ 3) e -^)|v;). 

It follows that 

£)(3) e -(«|«) = ^( 2 ) e -( 2 l u ). 
(The operator Z)£, 2 ' ) acts on the variable 2, while acts on the variable it.) [j 

The Laguerre functions are eigenfunctions of D^ 3 '': 

dW*£> = (2|m|+rM/)*^. 

(Actually the Laguerre functions are solutions of a system of differential equations and the 
above equation is one of them, see [Ricci-Vignati,1994].) 

In order to determine the operator Du acting on functions on C n such that T v = 
T V D^ we need some preliminaries. 

Lemma 5.5. 

n d 

tr(V^ B (x)) = J2( S J + ~ l))aj(s)<p s - £ .(x). 
i=i 

(In agreement with Lassalle's results [1998], p. 321, first line of (14.1).) 
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Proof. For t > we consider the following Laplace integral: 

f e -(x|v) e -*tr Vs(y ) A ( y )-f TO (^) = r n (s + p)^_ s (te + x). 

Taking the derivatives with respect to t for t = 0, one gets: 

- / e-W y hry<p s (y)A(y)-%m(dy) = Tc l (s + p) tr(Vp_ 8 (x)). 
By using Proposition 4.1: 

n 

tryy? s (y) = ^2aj(s)ip s+£j (y), 

3 = 1 

and since 

n „ n 

YV(s) / e- (x|y Vs+ £j (?/)A( ? /)-^m(rfy) = V oy(s)r n (s + + p)<p_ s _ e ,(x), 
,=i Jn ,=i 



one obtains 



tr(Vy-(x)) = - a 3 -(s) ^± g ± p) ^_ s - £j (x) 



-^«j( S )( s j - 4 fa" l))^-s- £ ,(^), 

i=i 



or 

tr(Vcp s {x)) = J2 a j(~ s )( s j + l( n ~ lj)<P*-ej(x). 

3 = 1 

In fact, by the explicit formula for Tn, 

Y n {s + p) = {2'K) N -\[Y{s J --{n-l)), 

3 = 1 

and 

r n (s + ej - + P ) _ r( 8j - + i- j(n-i)) _ d 

ra(s + p) - T( Sj -i(n-l)) ~ Sj 4 {U ^ U 

Lemma 5.6. 

d 

(Vp s (x),x 2 ) = J2( s 3 - |( r - 1 ))«j( S )¥ 7 s+ £j (^)- 
i=i 
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Proof. One uses the following formula: if F{x) = f(x 1 ), then 

(VF(x),x 2 ) = -(V/(x" 1 ),P(x- 1 )x 2 ) = -tr(V/(x -1 )). 
Taking f(x) = </?_ s (x), = <^ s (a;), one obtains: 
(V¥s(x),x 2 ) = -tr^-s^- 1 )) 

n ^ 
= - J3( _S J + 1 ))"j( S )V-s- £j (^" 1 ) 

n o! 
3=1 

Lemma 5.7. 

Dl 2) <Ps(z) = ^2(s j + u--(n-l))a j (a)<p 8+ej (z) 
3=1 

n <i 



For s = m — p, one obtains 

(D^^ m )(z) = J2{m J +u- -(J - l))«,(m - p)$ m+ej (z) 

3 = 1 
n ^ 

J=l 

(Theorem 2.2 in [0rsted- Zhang, 1995b] is an equivalent result.) 
Proof. By Lemmas 5.5, 5.6, and Proposition 5.1, we obtain 

D^ip s (z) = (z 2 ,V<p s (z)) - (e,V<p s (z)) + v(trz)<p s (z) 

n d n d 

= J2( S J ~^ n ~ 1 ))«j( s )v ? s+ £j (^) - + j(n - l))ctj(-B)(p a - E .(z) 

i=i j'=i 

n 

+ z/ $^«i( s )y 3 s+ ej (^) 

i=i 

J=l 3=1 
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Lemma 5.8. 

s£ 3 W«) 

n d d 

n 

+ l]«j( S )Vs+ ej (w). 

i=i 

Proof. Recall that the operator is defined by the relation C V D„ = D^C V . The 
modified Laplace transform of <p s is given by: 

r\rv f 
1 



n 



Tq(s + u + p)<f_ s _ v (z). 



We obtain, by using Lemma 5.7, 

£>£, 2 > (£^ s ) (z) = YTT^^n (s + v + p)x 



and this can be written 



+ XX S J +z/ ~ l))a J (s + z/)v?- s - ; ,- ej .(2)J, 



and 



- X> + 4< n - 1} )^ ( - S " y ) r n (s + y + p- ej -) ( ^- ej)( ^ 

E, d. , . r(s + p + v) . „ w . 

.J"' + " - 3 (» - i)W» + y ) rn (g +y +p+ ei ) (£ -^ )W - 



r n (s + i/ + p-ej) r( Sj - + i/-l- |(n-l)) ' 4 
Tn(s + p + z/) _ r( Si + i/-|(n-l)) 



Tnis + v + p + ej) r( Sj + v+l-$(n-l)) sj + u-^n-l)' 
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Finally one observes that ctj(s + v) = ctj(s). M 

End of the proof of Theorem 5.2 

(i) Let be the operator defined by the relation D„ T v = F^D^. We will show 
that D^f 1 = D v . The modified spherical Fourier transform T v is given by 

(•F„V)(s) = -p 7 ~~~ v , \ [ ¥>s(u)ij(u)A(u)%-%m(du). 

r^(s+ 2 +p) Jo. 

(3) 

Using the fact that the operator Dt, ' is symmetric with respect to the measure 
A(u) u ~~m(du) : one obtains 

D^\T^){s)=F„{D^){s) 

1 f Lp s -»(u)(DWiP)(u)A(uy-^m(du) 

(L> (3) 'ip s -»)(u)*jj(u)A( U y-™ m{du). 



r n (s+f +p) Jn 
1 



r n (s + f + p) J n 

From Lemma 5.8 it follows that 

D< 4) (JW)(s) = 



aj(-s) / y a - e .(u)%l)(u)&*~"(u)m(du) 

I n r N 

+ r 7 — 17— x / ^+ e »^(«)A5""(M)m((l«) 

J- n (s + 2 + Pj i=1 in 

= - - \ + 4 ( n - !)) + i - 1 - t {n ~ 1 aj(_s) 

(F v if>)(a-ej) 



j'=i 

rn(s + f - gj+p) 



™ ,r n ( s + § +sj + P ) 



As we saw in the proof of Lemma 5.8, 



rn(s + p+f) ^+1^-1-1(^-1)' 

r n ( s +f + p + £j ) ! d 

— j; . — ^ =Sj + 5i/- -(n - 1). 
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This shows that D^f 1 = D v . Since the Laguerre functions are eigenfunctions of the 
operator hence the Meixner-Pollaczek polynomials Qm are eigenfunctions of the 

operator D^f 1 = D v : 

D u Q^> = {2\ m \ + v)Qt ) . 



(ii) The identity in Theorem 3.1 can be written, with s = m + ^ — p, 

^ — n\&iJnA = A(p — iu\~ v (b~ ((p. — in} (p. 4- m) -1 ) = 



J2Qk } (m+ £ - p)$kH = A(e - w)-"Q m ((e - w)(e + w)- 1 ) = (-1)^2^ (w). 



(2) 

The right handside is an eigenfunction of D„ for the eigenvalue 2|m| + v. The action of 

(2) 

D v on $k has been computed (Lemma 5.7), hence 



(2|m| dmQtt (m + \ - p)* m (w) 

k 1 

= E d *Q± ( m + i - ( + ^ - ^' - x )K( k - p)**+<* h 

k j=l 
n , 



By equaling the coefficients of $>k(w) in both handsides, one gets statement (ii), observing 
that 2|m| + v = 2 Y^=i s j: since X^=i Pi = 0- D 

As in the one variable case, the operators (i = 1, 2, 3, 4) are essentially selfadjoint 
(with apporpriate domains). 

6. A Barnes type integral. - Let p be a positive measure on M with finite 
moments of all orders: for every m > 0, 



/ \t\ m p{dt) < oo. 
Jr 



For d > one defines the following probability measure on 

1 



(d) |V(xi, . . .,x„)| d //(^i) • ■■v(dx n ), 



with 



z i d) = / \V(x 1 ,...,x n )\ d p(dx 1 ) ...p(dx n ). 
A general problem is to evaluate its moments: 

m ^ = —(d) I x{xi,---,x n )\V{x u ...,x n )\ d p{dx 1 ) ...p{dx n ), 
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where x is a symmetric polynomial. Special attention have been made to the case 

1 f n 

and, more generally, 

F K (z 1 ,...,z K ) = —rjr / F[ \~\(z k -x j )\V(x 1 ,...,x n )\ d ii(dx 1 )...ii(dx n ). 
zl> Jr™ £S M 



k=l j = l 



For d = 2 the use of orthogonal polynomials leads to the evaluation of F\ (z) and 

Fk(zi, • • • 7 -^k)- Let (p m ) be the sequence of orthogonal polynomials with respect to the 

measure //, normalized by the condition 

Pm(t) = t m + lower order terms. 

Then 

Z< d > =n\h Q h 1 ...h n - U 

where 

h m = p m {t) 2 n(dt). 
Jr 

Proposition 6.1. — (i) For z eC, 

F\{z) =p n (z). 

(ii) For z-l, . . . , z K E C, 



F K (zi, . ..,z K ) = 



det (P"+j-i(*fc))i<j,fc<jr 

V(Z!,...,Z K ) 



(See [Szego,1975] p.27, formula (2.2.10) for (i), and [Brezin-Hikami,2000] for (ii). In 
point of fact (i) appeared already in [Heine, 1878].) 
If fj, is the measure on M given by 

/ f(t)fi(dt)= f 1 f(t)(i-t) a (i + tfdt, 

Jr J-i 
with a, j3 > — 1, then Zn is the Selberg integral and the evaluation of 
If 

® t (x) = -[a / x(xi,...,x n )\V(x 1 ,...,x n )\ d T[(l-x i ) a (l + x i ) (3 dx 1 ...dx n 
Zr V fj[ 

has been considered by Aomoto, Kaneko, and Mimachi. 
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Proposition 6.2. — (i) For z e C, 

F 1 (z) = Cp i J a ^ l3 \z), 

where pm ,f3 ^ is the Jacobi polynomial of degree m. 
(ii) For z-l,...,z k E C, 

F K ( Zl ,..., z K ) = CP{l^^\z u z K \ 



where Pm'^' 7 ' is the multivariate Heckman-Opdam Jacobi polynomial for the partition 
m = n K = (n, . . . , n). 

([Aomoto,1987] for (i), [Kaneko,1993], [Mimachi,2001] for (ii).) 

For d = 2, formulae (i) and (ii) agree with formulae (i) and (ii) in Proposition 6.1. 
In fact, for d = 2, there is a determinantal formula for the multivariate Heckman-Opdam 
Jacobi polynomials. 

Mimachi considered the following Barnes type integral, for k > 0, 

if ^ \ 

m ^ = -73) / x(Ai, , A») II l r (*^ + ?)l 2 n77TY^ dAl • • • dX ^ 

Zn ' JR n - =1 Z \ C \ tA )\ 

where x is a symmetric polynomial, and 



f n 1 

z { n d) = ni r (^+f)i 2 u7T^^---^ 



Proposition 6.3. — Define 

1 „ n n 1 

* W = ^> /„„ J> - **) II l r <^ + \ > l 2 RaF dA ' ' ■ ■ <"«• 

T/ien 

F 1 (,)=«(^) n ,3 K) (^), 

where qm is the one variable Meixner-Pollaczek polynomial of degree m. 

([Mimachi, 1999]. Actually this is a special case of Mimachi's result: <f> = A = | 
with Mimachi's notation.) 

In case d = 2, the reciprocal of the c-function is the Vandermonde polynomial: 



|c(zA)| 2 

27 



= 4v(\ 1 ,...,x n y. 



Then Proposition 6.3 is a special case of (i) in Proposition 6.1. 



We propose an alternative proof of Proposition 6.3. In our setting Mimachi's integral 
can be written 

9K(X) = 4 4) / X (A)|r^A + p+^)| 2 — l—dXL-.dXn, 
Jm n 1 \ c \ lA )\ 

with v = k + |(n — 1). Recall that the map 

7 :D(fi) ^P(C n ) s ", D^ lD 

is an isomorphism. Therefore, given x £ V(C n ) s ™, there is an invariant differential oper- 
ator D such that 

X(A) = lD (-i\- -), 

and, by Lemma 2.1, x is a l so the modified spherical Fourier transform of D^q (^q(u) = 
e -tr(u)y Therefore, by the Plancherel formula, 

m( X )=a^ [ (D$ )(«)*o(«)A(u)"-»m((iu). 
Jo, 

a) If x(A) = 7m(— ^A— |), then D = D m , the invariant differential operator such that 
a D {u,e) = $ m (w), and £)* («) = \E f m (-'u)e~ tr(u) . In that case 

Wt( X ) = (-l) |m| 4 3) / $m(u)e- 2tr ^A(uy-%du= (-l)l m l2-l m l(z/) m . 

b) We consider now the invariant differential operator D a (a G C): 

D a = A(u) 1+a A^A(u)- a . 
(See [Far aut-Koranyi, 1994] p. 294.) The polynomial 7d« has been computed (p. 296): 

7d«(s) = Y[(sj -a+ -{n- 1)), 

i=i 

and by Proposition XIV. 1.5, 

a, = £(-i)*("W>o<»-*> 

fc=0 ^ ' 
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((k) denotes the partition (1, . . . , 1, 0, . . . , 0) with 1 repeated k times). Therefore, by a), 



Mix) = (-l)^-^ (fy(a) {k) (v) {n _ k) 2 k . 



The coefficient (ct)^ can be written: 

(«)<*> = a(a - ^) (a - 2^) . . . (a - (fc - 1)^) = (^) [^a] fc . 

Hence 

fc=0 v 7 

By observing that 

± Qt+H.J' = [„]„£ g)(-i ) *(-„ )t ^_ I ^ 2 . 

k — k — 

as we saw in the proof of Proposition 1.1, we obtain the result. [} 
A natural question arises: is it possible to evaluate the integral 

f K ( Z1 ,..., z K )) = -75) / n n( z * - a ^ n + f )i 2 utxt^ ai • • • rfA - 

in terms of the multivariate Meixner-Pollaczek polynomials Q^^(Ci> • • • ? 00 ? For <i = 2, 
by (ii) in Proposition 6.1 and Theorem 4.3, it holds: 

F K (z 1 , ...,z K ) = CQ^ K) {iz u . . . , iz k ). 
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